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q ' Abstract 

O 
O . 

In this paper we study a class of inhomogeneous cosmological models which is 
^ . a modified version of what is usually called the Lemaitre-Tolman model. We 

Q | assume that we have a space with 2-dimensional locally homogeneous spacelike 

surfaces. In addition we assume they are compact. Classically we investigate 
both homogeneous and inhomogeneous spacetimes which this model describe. 
For instance one is a quotient of the AdS4 space which resembles the BTZ 
. black hole in AdS3. 

Due to the complexity of the model we indicate a simpler model which 
| can be quantized easily. This model still has the feature that it is in general 

inhomogeneous. How this model could describe a spontaneous creation of a 



o 



universe through a tunneling event is emphasized. 



o 
o 
o 

I. INTRODUCTION 

In the eighties several authors [|l]-|J suggested that the universe could have been sponta- 
neously created though a tunneling event. What the nature of such a tunneling event could 
be is still unsettled. Today's theories suggest however that shortly after the Planck era the 
universe was in an arbitrary state. A homogeneous and isotropic universe is believed to be a 
rather special configuration, so an arbitrary state suggests the universe was inhomogeneous. 
Recently the question of topology^ has been added into the discussion. What is the most 
probable topology of the universe? We know that there exist a lot more hyperbolic manifolds 
than elliptic and flat ones ||. Do we therefore live in a hyperbolic universe? 

In this paper we shall study a class of cosmological models which has 2-dimensional locally 
homogeneous compact spacelike surfaces. The 2-dimensional compact oriented surfaces can 
be classified in term of an integer g, the genus of the surface. There is one elliptic surface 
g = 0, one flat g = 1 and a infinite of hyperbolic surfaces g > 2. 
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We write the metric in a similar form to that of Kuchaf 
rators (KKTY) 



and Kenmoku and collabo- 



ds 2 



-N 2 dt 2 + Q 2 {dr + N r dtf + R 2 dtt 



(1) 



N, N r , Q and R are all functions of r and t, and dfl 2 is the metric the g-holed torusQ, T g . 
These metrics cover a great variety of different spatial topologies. In terms of the Thurston 
classification |TTJ, their covering spaces are for instance: E 3 , S 3 , H 3 , MxS 2 and R x H 2 . 
Thus a great variety of cosmological models are incorporated in the metric (|T|). Our main 
concern will however be the Lemaitre-Tolman models (or LT for short). In the LT models 
g = 0. This case includes the first three of the Thurston geometric models. The case where 
R' — Q' — N' — are the Kantowski-Sachs (g = 0) metric and their hyperbolic versions 
(g > 2) studied by Fugandes |I2| . Throughout this paper we will denote all the cases KxT 9 



as KS. This model will link the models studied in the paper by Nambu and Sasaki [I3| with 
for instance the homogeneous FRW models or the Einstein deSitter models. This model can 
therefore among others describe a quantum tunneling of a universe from a Schwarzschild 
Black Hole into a deSitter universe. 

We will choose the curvature of T g to be k = 1,0,-1 for g = 0, g — 1 and g > 2 
respectively. The Gauss-Bonnet Theorem then implies 



A g K = 4tt(1 - g) 

where A„ is the area of the g— holed torus T„. If we choose A\ = 47r we can write 



(2) 



Ag = A-KUg 



(3) 



where 



Tin 



1, 9 = 0,1 

(<?-!), 9>2 



(4) 



From the action of pure gravity^, Sq = J d^x(R — 2A), we obtain for these models 



S, 



G 



J dt J drn g - A^- 1 \ R(-Q + (QN r )')(~R + R'N r ) + ~Q(-R + R'N r ) 2 ^j 
+ N (-Q- l RR" - ^Q^R' 2 + Q- 2 RR'Q' + ^Q(n - AR 2 ] 



(5) 



2 Note a slightly different notation compared to the notation used by KKTY and Kuchaf 

3 Or more correctly: dQ, 2 is a metric on the (/-holed torus, T g . Since the Teichmiiller space for 
g > 1 have dimension greater than 1, there are an uncountable number of non- isometric metrics 
on Tg, g>l. 

4 In this paper we have used the convention: c = Gn = h = ks = 1- 
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where the prime and dot denote the derivative with respect to r and t respectively. A 
cosmological constant is also included. We will also consider a matter term or source term 
of the form: 



Sd = —4:im g dt drNp 



(6) 



where p an independent function. We introduce this source function because it is one of the 
simplest ways to generate inhomogeneities. In the presence of classical dust this function p 
could be related to the dust density p by the equation p = QR 2 p. Thus in this description p 
is not a dynamical matter field, we have basically solved the matter equations and reinserted 
the solutions into the LagrangianQ . 

The total action is now assumed to be S = Sq + So- The canonical momenta are: 



P< 



Q 



P 



R 



n, 



[R(Q-(QN r )') + Q(R-R'N r )] 



N 



Through a Legendre transform we can bring the action into a canonical form: 



S = J dt J dr[P R R + P Q Q - (NSj + N r f, r )} 



where 



1 1 Q 2 

n g R 2n g R 2 v 



+ n 



RR" 

S) r =R'P R - QP'r 



RR'Q' 1 R 



r2 



1 -kQ + ^AQR 2 + 4vrp 



Q 



(7) 
(8) 



(9) 
(10) 



II. THE CONSTRAINTS IN THE KSLT MODEL 

From the canonical form of the action we can readily write down the energy constraint 
and the momentum constraint: 



Sy = 0, S) r = 

It is also useful to study a special combination of these two constraints, 



R Pq 

ft —fir 

Q RQ 



n g (M' - m) = 



(11) 



If we want to treat dust as a dynamical field we get a term in the Lagrangian that is linear in 
the conjugated momentum of the dust field. See for instance [ p| , |T5| ]. 
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where M is called the mass function and is defined by 



M 



P 2 



Rn 2 g 



and m is defined by 



m 



RRT 



drA-rrp 



kR + ^_R 3 



R!_ 
Q 



Inserting the relation p = pR 2 Q we see that in the case k = 1 we can interpret the function 
m as the mass of the dust inside a spherical shell of radius R. 
If we calculate the equal-time Poisson algebra we obtain: 



{$)r(r),Sj r (r*)} = f> r (r)6\r - r*) - f) r (r*)8'(r* - r) 

{fi{r),S) r {r*)} = &(r)6(r - r*) + Sj{r)8'(r - r*) 
{S)r(r),S) r (r*)} = ftr(r)5'(r - r*) - ?> r (r*)8'(r* - r) 

{i5(r),^(r*)} = Qtr)- 2 ^)^ - r*) - Q(r*)~ 2 fi r (r*)8'(r* - r) 
{M{r),fi r (r*)} = M'(r)5{r - r*) 

{M(r), ^(r*)} = -^(r)- 3 ^'^)^^)^^ - r*) 
{M(r),M(r*)} = 



(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



Thus according to the above equations the mass function is a constant of motion in a weak 
sense. 



III. CLASSICAL SOLUTIONS 

For simplicity, let us work in the gauge iV r = 0. But for the time being N is still arbitrary. 
The constraint $) r can easily be solved in quadrature. Using the expressions for the canonical 
momenta, eq. and||, combined with the constraint equation 10, and simplifying yields the 
differential equation: 

• N' 

R'Q - QR' + RQ— = 
which has the solution R' = FQ. The function F satisfies 



• N' 

R'F - FR— 

N 







For the present time we have to be a bit careful when we solve this equation. If R! = 
which happens in the KS case, the momentum constraint yields N' = 0. For R' ^ we can 
solve the equation for F: 



F { r ,t) = f{r)exp 




RN' . 

dt 

R> N 
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and f(r) is an arbitrary function of r. We can add these results together and keep the above 
expression for F(r,t) but with the additional restriction f(r) = whenever R! = 0. We will 
ignore the possibility that the term ^ in the exponential function causing trouble^]. Thus 
if there exists a t* and a r* such that R'(t*,r*) = 0, then F(t,r*) = for all t. 
The mass function is always well defined and turns into: 

M = ^ + V K - F2)R ~^ = m (19) 

where m(r) = f drAnpF. If F = the mass equation ([II]) is separable and can be solved 
in quadrature. 

A. Time independent solutions 

Let us assume R = 0. This implies F = 0. Solving the mass equation with respect to F 
yields: 

„ 2m(r) A o 

F = k ^ ^ 2 

R 3 

In order to find the lapse iV we have to use the equation for Pr\ 

-^ P Q = ^{PqJ drN *} (2°) 

1 Po (\ R' 2 I, . n2 , N'RR'\ 

+ % o7^-o(^- A ^) + ^77^ (21) 



2n g R 2 y \2Q 2 2 y ' ' NQ 
Inserting F 2 into the equation for Pq yields the following equation for the lapse iV: 



N' _ 1 / 1 k — AR 2 

~N~2\R + kR — 2m — #i? 3 



o + — - a^ K (22) 



As we see, we might as well (at least locally) consider R as our new radial variable, i.e. 
R = r. If we do so, the solution to the above equation can be written as: 

r 



where 



dr 



k — Ar 2 



kt — 2m(r) — |-r 3 
and rj{t) is an arbitrary function of t. Thus 



6 We get a "^"-expression for £L. 
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ds 2 



-r, 2 (t)-e I dt 2 + 



K — 



dr2 + r 2 dn 2 

2m(r) A 2 9 



(23) 



- 3 r 



wherever k — 2ml - r ' — 4r 2 > 



If for instance we have constant dust density, i.e.[] m(r) = \i + |-yrpr 3 the integral may 
be evaluated in terms of elementary functions [TH . 



Spatial topologies in AdS^ space 
For all the deSitter spaces (m = 0) the integral I is simply: 

T ^ ( A 3 

i = In \ kt r 

V 3 

for physical spacetimes. We notice that k < is not compatible with a A > 0. But 
for A < we may have any k. All of these spaces will have spatial sections which have 
constant negative curvature, i.e. they belong to the Thurston classification H 3 . The cases 
k = 0,-1 are not compatible with compact spatial sections, but k = 1 may have compact 
spatial sections if we demand that the metric |23] is defined locally in r. If we demand that 
the metric should be globally defined then we obtain the space H 3 itself. Unfortuneatly a 
classification of the possible compact quotients of H 3 is still far from being established. Let 
us for simplicity's sake set A = —3. 

• k = 0: This is the non-compact multiply connected space investigated by Sokoloff and 
Starobinskii jl^]. This space is best illustrated in the Poincare half space modelQ. If 
we do the coordinate transformation r = — - the metric turns into 

z 

ds 2 = ir(-dt 2 + dx 2 + dy 2 + dz 2 ) 
z l 

which is manifest on Poincare half space formQ. We identify points along two pairs of 
planes parallel to the z-axis. It is easily seen that this space does not allow compacti- 
fication in the z (or r) direction. 

• k = — 1: This section of the hyperbolic space can be illustrated as follows. Draw two 
4g-gons on the boundary of the Poincare ball model. These polygons can be placed so 
that they are symmetric to each other through the origin. Connecting the polygons 
using geodetic planes we obtain a space that look like an apple core(see figure H). Also 




7 Note that in general p ^ p. We define p so that it is not a funtion of t, while p in general could 
be a function of t. 

8 See figure [j]. 

9 Doing a Wick rotation of the time component, r = it, we obtain the "Euclidean" section of AdS4 
space. Here we see that the Euclidean section is the 4-dimensional hyperbolic space H 4 . 
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this case will not allow a compactification in the r direction. Actually we have cheated 
a little bit, because the space contains a horizon with topology as a g-torus, much like 
that of a black hole 1171 . The of this horizon is 



127T 

9 ~ JK\ Ug 

Thus apparently we can relate to it the Bekenstein-Hawking entropy: 

S BH = - A J T dA = -n 9 



This metric is very similar to the BTZ black hole in AdS3 [flgl - Including a mass 
term in this case is straight forward. Note also that in the above case the "Black 
Hole" entropy is a feature of the negatively curved compact spaces T g . To obtain 
the AdS4 apple we have actually taken two such spaces and glued them along their 
horizons to obtain a topological complete space. Note also that the horizon is the only 
r = constant surface which is totally geodesic. To ensure regularity of the total space 
after we have performed a gluing we have to glue along totally geodesic surfaces. 

If we had demanded that all of these cases should be simply connected, they would have 
yielded the same space H 3 . But since we have assumed that our locally homogeneous 2- 
dimensional surfaces are compact all of the above cases yield topologically different spaces. 
Again |2(J we see that assuming compact 2-dimensional locally homogeneous spaces ensure 
that the different cases will yield different solutions. 



B. Kantowski-Sachs-like solutions 



In the Kantowski-Sachs case we have B! — 0. This implies F — 0. The mass equation 
will now yield: 



RR 2 1 A 



(25) 



where /J, is a constant. Inserting this into the energy constraint we obtain an equation similar 



to the equation |22 



9. 

Q 




SitpR 2 - k 



±R 3 + 2^-kR 



R 



(26) 



Again we see that we can equally well consider R to be R = t. This equation is very similar 



to the equation |22|. This is not surprising since if we smoothly continue the Schwarzschild 
metric to the region inside the black hole horizon we obtain the KS metric inside the horizon 
T3| . Smooth continuing the deSitter solution to the region outside the deSitter horizon gives 
the KS metric |21|. We also note that the t <-> r symmetry is manifestly broken with the 
inclusion of classical dust. Due to the conservation of the rest mass of the dust, we would 
obtain J s d 3 xpQR 2 = V C M. where V c = J E d 3 x is the coordinate volume if the spatial sections 
S are compact. 
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We re-write these equations for a general dust distribution. We assume that the topology 
in the r-direction is that of a circle S 1 . If < r < £ then we can write 



n=l 



a n (t) cos 



2imr 



+ a_ n {t) sin 



2irnr 



(27) 



We use the same trick to write (|t 2 + ^ — k) 2 Q = G where 



G(r,t) = J2 



n=0 



( 2nnr\ 
c n (t) cos I — - — I + c- n (t) sin 



/ 2imr\ 



(28) 



Since the cosines and the sines form a complete and orthogonal set of functions on the 
circle the equations for c n , n^O are 



c n (t) = 



47ra n (t) 



t(¥ 2 + f -«) 



(29) 



The equation for c (t) is 



c (t) 



A-kM 



(30) 



can be integrated at once by elementary calculus. 

In the absence of dust p = this reduces to the ordinary inner/outer solutions of the 
Schwarzschild-deSitter spacetime. This is just the homogeneous solution of the KS metric. 
If there are inhomogeneities we have to solve the equations for c n (t), n^O. Unfortunately 
we then have to know the evolution of the matter content. If we on the other hand knew 
the metric, i.e. the functions c n (t) we could uniquely determine the matter distribution. 



IV. LT SOLUTIONS 

So far we have only investigated time-independent and KS-like solutions. We saw that 
they represented different sections of a more general spacetime. Let us now look at the LT 
solutions, which we can obtain if we choose the universal time gauge N — 1. The coordinates 
now follow the collapsing/expanding dust. 

If we choose N — 1, F — F(r) will follow. We can now write the mass equation as an 
"energy" equation: 

±R 2 + V(r,R) = E(r) 
where the "potential" V and the "energy" E are given by 

R 6 
E = -\{k-F*) 
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This energy equation may be integrated and solved exactly. A summary of the results 



and a qualitative description of the physical meaning of the solutions is given in |22]. The 



solutions may be written in terms of the Weierstrass' elliptic functions [p3 [. The actual 
expressions are not very informative unless the reader has massive knowledge of these elliptic 
functions. However a lot of qualitative information can be extracted from simple classical 
considerations. 

We will however solve the energy equation in the big and small R limit. 

Large R: Neglecting the m term 

If we neglect the m term in the energy equation, we get a simple and easily solvable 
equation. The equation becomes: 

R 2 - H 2 R 2 = 2E 

where H 2 = ^ (assuming A > 0) This equation can be solved in quadrature, and the 
solutions are deSitter-like solutions^ 

^P^sinh(tf(t-to(r))), E >0 
R(r,t) = {e H <t- t °W, E = (31) 

^ 2 g (r)l cosh(ff(t-t (r)), E<0 



Small R: Neglecting the A term 

If the A term is neglected, we end up with a energy equation which looks like: 

•9 2m „ 

R 2 = 2E 

R 

Introducing a parameter G we can also solve this equation in quadrature: 




m 

2E 



(cosh 0-1) 
(sinhG-6) f>£ >0 ( 32 ) 



(2E)i 

R= (Jm) (t - t (r))l , E = (33) 




2 

3^(1- cos G) 
^^(G-sinQ) 




(34) 



From these solutions we notice that all the small R solutions goes as 

1 2 

R ~ (^ m ) 3 (t — to(r))3 as R — ► 0. In the special case to(r) = constant and F 2 = 1 — kr 



2 



the equations |32], |33] and |34] reduces to the well known matter dominated homogeneous FRW 
solutions. 



10 For a rich and qualitative description of various deSitter models see for instance [24] 
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The behavior of a general solution 

The classical solutions will move on level curves of the energy function E(R, R) = ^R 2 — 
^ — ^R 2 , since the total energy E is independent of t. In figure [| the level curves of a typical 
energy function are drawn. If A > there will exist a saddle-point of the energy function. 

T his saddle point will be at R = (^p) 3 , R — where the energy function will have the 
value E s = — i(9m 2 A)3. The saddle point solution is static, and is as a matter of fact the 
Einstein static universe^. If E < E s , the solutions fall into two distinct classes: 

1. Schwarzschild-like solutions: These solutions expands, but they do not possess enough 
energy to escape the gravitational collapse, so they end as black holes. If m is constant 
these solutions are those of a Schwarzschild black hole in Lemaitre coordinates. 

2. deSitter-like solutions: solutions where the universe evolves approximately as that of 
deSitter solutions with positively curved hypersurfaces. 

If E > E s the (test) matter has enough energy to escape the gravitational collapse (expand- 
ing solutions) or enough energy to prevent the gravitational repulsion from the cosmological 
constant (contracting solutions). 

NOTE: 

All of these solutions are also valid for the case F = the KS case. But to obtain the 
KS solutions it is not enough to solve just the mass equation. 

V. THE QUANTUM MECHANICAL LT MODEL 

From here on we will only consider the case where the cosmological constant is non- 
negative. As we saw from the classical case, if the "energy" is small enough the solutions 
are confined to be either Schwarzschild-like or deSitter-like. These two regions are separated 
by a classically forbidden region. From classical quantum mechanics we have learned that 
this does not have to be the case for a quantum system. A quantum mechanical system 
may tunnel though classical barriers. In our case this means that a small Schwarzscild-like 
universe may tunnel though the potential barrier to become a deSitter-like universe! This 
has striking consequences. If we imagine that a "mother" universe could spontaneously 
create Schwarzschild-like black holes of less than a Planck-length, there will be a non-zero 
probability that one of these will tunnel through the potential barrier, and end up in an 
ever-expanding deSitter state. These considerations are of course rough and based on clas- 
sical quantum mechanics. Whether this scenario is probable in the yet unknown theory of 



This is easily seen if we define the new radial variable to be R. This can be done since we 
have to assume that m'(r) > on physical grounds. In the case m'(r) = 0, the metric become 
degenerate. 
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Quantum Gravity remains to be proven. However we will use a simplified version of the WD 
equation and estimate the form of the wave function in various limits. We will from now on 
consider only the spherically symmetric case, i.e. k = 1. 



A. The canonical Dirac-quantization 

The Dirac quantization procedure is to let the constraint equations become operator 
equations of the wave function: 

m = (35) 

S)r^ = (36) 

The first of these equations is the Wheeler-DeWitt equation. In addition we have the 
momentum constraint equation. These two equations are functional differential equations 
and thus in general very difficult to solve. KKTY managed in the case p = to find a couple 
of special cases where the wave function was exactly derived and interpreted [53] . Since the 
mass equation is a special combination of the momentum and energy constraint, the mass 
equation itself can be considered as a constraint. The wave function has to obey: 

M# = m^I (37) 

The wavefunction has to be an eigenfunctional of the mass-operator M with eigenvalue 
m(r), and from the classical interpretation we have to assume that the eigenvalue m(r) has 
to obey m'{r) > 0. 

In the WKB approximation we can integrate the Hamilton- Jacobi equation if we choose 
a particular and simple path ||26|| . We will not use these results. Neither will we give any 
attempt to solve all three constraint equations simultaneously in the "correct" quantum- 
mechanical manner. To simplify the equations we will do the following: 

• Solve the momentum constraint classically. When we solve it we will assume that 
B! 7^ 0, thus the KS case is excluded. Solving the momentum constraint will reduce the 
number of canonical variables, so we end up with only R and its canonical conjugated 
momentum n^. 

• Choose the gauge N' = so that we will not have any problems with the function F 
being dependent of t. 

• Find the new Hamiltonian which reproduces the classical equations. 
The reduced action which reproduces the classical equations is: 

S= [ dt [ dr— --^- + -R(l - F 2 ) - -R 3 - m(r) 
J J F 2 N 2 2 J 6 K ' 

Calculating the canonical momentum 

5S ( N\ ' RR 



SR \F) N 2 
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doing a Legrendre transform to bring the action onto a canonical form, and removing bound- 
ary terms, we arrive at the Hamiltonian which reproduces the classical equations: 



H 



i F 4 n| i 



F 2 



+ -R 6 + m r 
o 



(38) 



The constraint equation will now yield the correct classical equation; the mass equation. 
Calculating the equation for IIr, we find the identity, M = 0, which tells us that the mass 
function is independent of time. 

The wave function has to satisfy the constraint equation: 



d_ 

Or 



m = o 



(39) 



where 



k 1 Fi 
£ = - 



2 RF' 2 5R 2 



— R 3 + m(r) 
6 



A possible constant, due to the derivative with respect to r is assumed to be absorbed 
into the function m. 

This simplified version of the LT model is not only simpler but the action has the great 
advantage that it reduces to the deSitter models studied by Vilenkin and Hartle-Hawking 
if we demand that m = 0, R = a(t)r and F 2 = 1 — /cr 2 .0 The solution most likely 
to describe the tunneling behaviour emphasized in the beginning of this section, is the 
tunneling wavefunction suggested by Vilenkin. 



B. The tunneling wave function in the WKB approximation 

In the WKB approximation we assume that the wave function has the form ^wkb 
exp(±zS'), to the lowest order we get the Hamilton- Jacobi equation: 



5RJ 



F' 2 



2mR-R 2 (l 



F 2 ) + ±R* 







(40) 



If we assume that S = J a(r)dr the resulting equation will be the Hamilton- Jacobi equation 
for a point particle with action <r(r) (r is only a parameter). In the Hamilton- Jacobi equation 
the functional 5* turn out to be the action at the classical level. Since the classical action can 
be written as an integral over r the assumption S = j a(r)dr is therefore reasonable at the 
lowest order WKB level. This separation of the action is also a feature of the dust, in the 
sense that dust does not self- interact. Because the matter does not self-interact the matter 



12 Inserting m = |-7r/5r 3 , R = a(t)r and F 2 = 1 — kr 2 where p and k are constants into the action, 
the r dependent part can be factored out, and the integral may be performed. The actual integral 
obtained is: f dr , r , ,, = in the case of k > 0. 
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equations do not depend on neighbouring points. We can interpret the action a as the ac- 
tion of a point particle moving in a potential V(R) = ^ [—mR + |(1 — F 2 )R 2 — f-R 4 ] 
with zero energy. The WKB wavefunction ipwKB for the point particle can then be 
written i{jwkb = exp(±ia). The two WKB wavefunctions can therefore be related by 
^wkb = exp(J drlnipwKB)- Finding first the wavefunction if) we can then relate its WKB 
approximation with ^wkb through ^wkb = exp( f dr In iPwkb) ■ 

We imagine that the universe point starts off at R = and the behavior of the wave 
function of such a particle is dependent on the form of the potential. The zeros of the 
potential indicate classical turning points, separating classically allowed/forbidden regions. 
The potential V(R) has one zero at R = and possibly two more for positive values of R. 
Whether there are two, one or zero, depends on the entity: 

9m 2 A 
(1-F 2 ) 3 

Iff < I < 1 then there will be two positive roots, r\ and r 2 . The region r\ < R < r 2 
(assuming r\ < r 2 ) is a classical forbidden region, thus the wavefunction is exponential in 
some way. Iff Z = 1 these two roots have emerged into one, and the case I > 1 and / < has 
no positive roots; the classical particle rolls down the slope to infinity. 

One case which is very instructive to look at is the case where F 2 = 1 — kr 2 and the dust 
density is constant: m = ^npr 3 . The entity I will then be a constant with respect to r: 

? _ 16n 2 p 2 A 
k 3 

We can interpret this quantum mechanically as if we imagine a spontaneously created uni- 
verse "bubble" with high curvature, and a relatively small cosmological constant and dust 
density. It is most likely to vanish before it has the opportunity to enter a deSitter phase. 
Let us look at this from another point of view. Imagine that we can only observe "big" uni- 
verses, universes at scales much larger than the Planck length, or essentially just universes 
that are in the expanding deSitter stage. If we could in some way sit outside and see the 
flow of universes in expanding deSitter stages, we would find that most of the universes with 
low A and p also have low hypersurface curvature. 

Let these thought experiments be for the moment because we do not actually know how 
to interpret the wave function of the universe or how to calculate it. Let us however try to 
estimate the behavior of the wavefunction is this model. 

The most interesting case seems to be < I < 1, so let us just look at this case. There 
are essentially three different regions with different behavior. 

1. Region I: Small R. We can neglect the i? 4 -term. 

2. Region II: Intermediate region. This is dominated by the part where the wave function 
has an exponential behavior. 

3. Region III: Large R. We can neglect the i?-term. 
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Region I 

We neglect the i? 4 -term in eq. [H] and complete the square with a new variable —^=A 



(1-F 2 UR — — g- . Then the Schrddin ger equation for ib reduces to a Harmonic oscillator- 

V ' (1-F 2 )3 

like equation: 



+ (2p + 1 - A 2 )</> = 



where p = h ( — \Hhn — — 1 | . The solution of this equation can be written in terms of the 
Parabolic cylinder functions D p {z): 



'd p (±V2\) 
D_ p _ 1 (±iv^A) 



In the case of p = n an integer, the D p (z) can be written in terms of the Hermite polynomials 
(as one would expect): 

D n {z) = -2-2e~H n (—) 
The index p can be any number, in contrast to the 1-dimensional Harmonic oscillator. 



Region III 



We neglect the i?-term in eq. 4C. This is equivalent to setting the mass function m 
equal to zero. We are then left with a deSitter space-time. The wavefunction for such 
space-times has been calculated using the WKB approximation. Most interesting to us is 

the tunneling wave function of Vilenkin. Through the rescaling: r] = 



1 - F 2 )*R we 

can bring the equation to the same form as that of Vilenkin [E7I] , with a scalar field potential 



V 



F 2 A 



3|F'|(1-F 2 )2 



We use this result to obtain: 



1> 




R 2 A 
3(1-F 2 ) 



- R 2 A 
-3(1-F 2 ) 



R 2 A 
3(1-F 2 ) 

R 2 A 
' 3(1-F 2 ) 



< 1 
> 1 



(41) 



Let us assume that m(r) = and F 2 = 1 — kr 2 where k is assumed to be positive. We will 
now use the tunneling solution eq. HTI to find out what is most likely, the creation of a highly 



curved or a less curved space. If we calculate the ratio %*^\^^ where R 2 >> -|(1 — F 



get: 



we 



P k 



\[/*v]/ 



R=R 
R=0 



exp 



A 



d(sm9)G(8) 



(42) 
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where G is a real function independent of k which satisfies G{6) > 0,W6. We thus have 
to conclude that to lowest order in the WKB approximation it is equally probable for a 
space with small spatial curvature to tunnel through the classical barrier than for a highly 
curved space. Comparing with our previous statement, we se that larger curvarture itself 
increases the probability amplitude. But smaller curvature increases the volume by exactly 
the amount so that it cancels the contribution from the curvature. 



Region II 

Here the wave function will be of exponential behavior. But the actual shape is strongly 
dependent on the boundary condition. To find the wave function in this region to lowest 
order, the best is perhaps to use the function from Region I and III and cut and glue them 
together. 



Connection with other work 



Nordbury p8| , |29|| studied the FRW model with different types of perfect fluids. The 



tunneling picture which we have described here seems to agree with his results for a closed 
universe (1 — F 2 > 0). However in the case of open (k = —1) FRW-universes he concludes 
that these models can not have been created through the tunneling picture. From our 
considerations however, this can not be ruled out, on the contrary, open universes do seem 
to tunnel from a R = state a lot easier than closed ones. 

Also an article by Atkatz and Pagels || seems to indicate that open universes can not 
be created by a tunneling event. A key argument used in both of these articles is that 
the transition amplitude is suppressed by the infinite integral over the three-dimensional 
hypersurfaces. However, a more recent article by Coule and Martin |30| used the fact that 



there exists k — — 1 compact spacelike hypersurfaces. Since they are compact they will 
have finite volume. Their conclusion is opposite of that of the previous two articles: The 
creation of a open k — — 1 FRW universe is more likely than a closed k = 1 universe through 
a tunneling event. Our case is however more subtle, since the hypersurfaces are no longer 
homogeneous. 



C. Inhomogeneity Vs. Homogeneity: Tunneling Probabilities 

As we have investigated the WKB wavefunctions in this model, we might wonder whether 
an inhomogeneous universe is more likely to be created through a tunneling event than a 
homogeneous one. As an inhomogeneity measure in this model it is very useful to use the 
Weyl curvature tensor C Q/ g 7 5. All of the FRW models are conformally flat, thus for the all 
of the FRW model the Weyl tensor will vanish. In our case the opposite is also true, if 
C a p 1 s = and N' = in the LT model then the spacetime is spatially homogeneous. If we 
define a mean dust density function p by: 

4 

m (r) = -npR 3 (43) 
3 
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we can write the Weyl curvature invariant in the LT model as: 

16 2 

C a ^ 5 C a ^ 5 = -^n 2 (p-p) 2 (44) 

Thus this space-time is homogeneous if and only if p = p. 

The inhomogeneity in this model is therefore almost completely characterized by the 
mass function m(r). The homogeneous mass function will go as m^(r) = Kr 3 , where K is 
a constant. Investigating only closed universes, we fix the constant K so that m,h(r max ) = 
m(r max ). The most physical interesting universes that these LT models describe will have a 
larger amount of dust consentrated in the inner regions than the outher regions. The most 
extreme inhomogeneous mass function is the Schwarzschild mass function, ms c h{f) — In- 
constant. Thus it is very reasonable to assume that we have a mass distribution that obeys: 

m Sc h(r) > m(r) > m h (r) (45) 

If we compare the WKB tunneling potential V(R) = [—mR + |(1 — F 2 )R 2 — f-R 4 ] for 
these different mass distributions we see that we get 

Vsch(R) < V(R) < V h (R) (46) 

keeping all other parameters fixed. Thus the more inhomogeneous the space is the smaller 
the tunneling barrier will be. In the WKB approximation this means that the probability 
of the creation of an inhomogeneous universe is larger than for a homogeneous universe. 
More specific the spontaneous creation of a Schwarzschild-deSitter like universe appears to 
be more probable than a homogeneous FRW universe. We also have be aware of that the 
matter Lagrangian describes dust. Dust is an idealised form of matter that has no internal 
pressure. As the dust density becomes high, either in an early stage of the universe or in 
large local inhomogeneities this assumption of a pressure-free matter Langrangian becomes 
somewhat artificial and unphysical. In these situations we know that the internal pressure 
becomes very important and even quantum effects from the matter fields are essential in the 
description. In this sence we may say that these calculations more or less are "valid" only 
with small local inhomogeneities or they describe the effect from Gravity itself to produce 
inhomogeneities. This point is also emphasized by Conradi pl,p2|. Conradi discusses the 



wave function of an anisotropic KS space with a dust Lagrangian. He argues that since 
dust interacts only locally, in the sence that it does not self-interact, this local "tunneling" 
process cannot be interpreted as a tunneling effect. Our results should therefore perhaps be 
seen from a more phenomological point of view. The matter Lagrangian used in this paper 
is one of the most simple source field that give rise to inhomogeneities. The source field used 
is also interesting because it connects the Schwarzschild spacetime with the FRW spacetime 
in a smooth way. 

VI. CONCLUSION AND SUMMARY 

In a Hamiltonian description of General Relativity there has been a consensus among 
physicists and cosmologists that the issue of topology is essential and necessary in the de- 
scription ]33| , 34}] . This has also been emphasized in this paper. At the classical level we saw 
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that apparently different cases would lead to the same spacetime if the symmetry spaces 
were infinite and simply connected. Assuming compactness and multiply connectedness this 
pathology was removed and yielded some interesting spacetimes. 

In the Quantum mechanical case we provided a simpler model which was in general in- 
homogeneous, but removed many of the algebraic difficulties in the full KSLT model. It was 
emphasized how this model could describe a quantum creation of a inhomogeneous universe 
through quantum tunneling. We also showed that under some reasonable assumptions the 
creation of an inhomogeneous universe was more probable than a homogeneous FRW uni- 
verse. In the absence of interacting matter this is very reasonable because there are a lot 
more inhomogeneous configurations than homogeneous ones. 

Even though the full theory of Quantum Gravity is far from being established, it seems 
like Quantum Cosmology could describe many of the features of the initial universe. In 
the recent decade the interest for topology in cosmology has expanded the application for 
Quantum Cosmology. Even in classical gravity the inclusion of non-trivial topologies has 
yielded new and surprising results. 
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TABLES 



A 


K = 1 


K = 


K, = — 1 


+ 


KS— S-dS— KS 


KS 


KS 





KS— Sch 


KS 


KS 




KS— AdS 


KS— AdS 


KS— AdS 



TABLE I. Summary of the classical solutions. The ones to the left are the interior solutions, 
to the right are exterior solutions. We also use the term KS in all cases where the spatial topology 
can be written as S 1 xT g . Other abbreviations: S-dS: The Schwarzschild-deSitter family, Sch: The 
Schwarzschild solution, AdS: Anti-deSitter (or quotients of). 
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FIGURES 




FIG. 2. The AdS4 apple. This figure is not drawn to scale. 
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